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General relation is derived which expresses the fidelity of 
quantum dynamics, measuring the stability of time evolution 
to small static variation in the hamiltonian, in terms of ergod- 
icity of an observable generating the perturbation as defined 
by its time correlation function. Fidelity for ergodic dynam- 
ics is predicted to decay exponentially on time-scale oc 8~ 2 , 
8 ~ strength of perturbation, whereas faster, typically gaus- 
sian decay on shorter time scale oc <5~ x is predicted for in- 
tegrable, or generally non-ergodic dynamics. This surprising 
result is demonstrated in quantum Ising spin- 1/2 chain peri- 
odically kicked with a tilted magnetic field where we find finite 
parameter-space regions of non-ergodic and non-integrable 
motion in thermodynamic limit. 

PACS number: 05.45.-a, 03.65.Yz, 75.10.Jm 



The quantum signatures of various types of classical 
motion, ranging from integrable to ergodic, mixing and 
chaotic, are still lively debated issues (see e.g. Qj). Most 
controversial is the absence of exponential sensitivity to 
variation of initial condition in quantum mechanics which 
prevents direct definition of quantum chaos ||. However, 
there is an alternative concept which can be used in clas- 
sical as well as in quantum mechanics || : One can study 
stability of motion with respect to small variation in the 
Hamiltonian. Clearly, in classical mechanics this concept, 
when applied to individual trajectories, is equivalent to 
sensitivity to initial conditions. Integrable systems with 
regular orbits are stable against small variation in the 
hamiltonian (the statement of KAM theorem) , wheres for 
chaotic orbits varying the hamiltonian has similar effect 
as varying the initial condition: exponential divergence 
of two orbits for two nearby chaotic hamiltonians. 

The quantity of the central interest here is the fidelity 
of quantum motion. Consider a unitary operator U be- 
ing either (i) a short-time propagator, or (ii) a Floquet 
map U = Texp(— i f£ drH{T)/K) of (periodically time- 
dependent) Hamiltonian H (H(t + p) = H(t)), or (iii) 
a quantum Poincare map. The influence of a small per- 
turbation to the unitary evolution, which is generated by 
a hermitean operator A, Us — U exp(— iAS), S being a 
small parameter, is described by the overlap (ips(t)\ip{t)) 
measuring the Hilbert space distance between exact and 
perturbed time evolution from the same initial pure state 
\ip(t)) = ?7'|V')j \*Ps(t)) — UsW)i where integer t is a dis- 
crete time (in units of period p) Q. This defines the 
fidelity 



F(t) 



(1) 



where the average is performed either over a fixed pure 



state (.) = (ip\.\ip), or, if convenient, as a uniform average 
over all possible initial states (.) = (1/TV) tr (.), TV being 
the Hilbert space dimension. The quantity F(t) has al- 
ready raised considerable interest, though under different 
names and interpretations: First, it has been proposed 
by Peres |3| as a measure of stability of quantum motion. 
Second, it is the Loschmidt echo measuring the dynamical 
irreversibility of quantum phases, used e.g. in spin-echo 
experiments || where one is interested in the overlap 
between the initial state \tp) and a state U7 f7*|^) ob- 
tained by composing forward time evolution, imperfect 
time inversion with a residual interaction described by 
the operator Ad, and backward time evolution. Third, 
the fidelity has become a standard measure characteriz- 
ing the loss of phase coherence in quantum computation 
Q . Fourth, it was used to characterize "hypersensitivity 
to perturbation" in related studies Q] , though in different 
contexts of stochastically time-dependent perturbation. 

The main result of this paper is a relation of the fi- 
delity to ergodic properties of quantum dynamics, more 
precisely to the time autocorrelation function of the gen- 
erator of the perturbation A. Quantum dynamics of fi- 
nite and bound systems has always a discrete spectrum 
since the effective Hilbert space dimension TV is finite, 
hence it is non-ergodic and non-mixing [|]J^]: time cor- 
relation functions have fluctuating tails of order ~ 1/TV. 
In order to reach genuine complexity of quantum mo- 
tion with possibly continuous spectrum one has to en- 
force TV — > oo by considering one of the following two 
limits: quasi-classical limit of effective Planck's constant 
H — ► 0, or thermodynamic limit (TL) of number of par- 
ticles, or size L — > oo. Our result is surprising in the 
sense that it predicts the average fidelity to exhibit ex- 
ponential decay on a time scale oc 8~ 2 for ergodic systems 
(i.e. such that the integrated time auto-correlation of A 
is finite), but much faster, typically gaussian decay on 
a shorter time scale oc S^ 1 for integrable and general 
non-ergodic systems (i.e. such that time averaged auto- 
correlation of A is non- vanishing) . Our theory on fidelity 
is very general and can be extended to any perturbed uni- 
tary evolution, either in quantum, quasi-classical, or even 
classical (Liouvillian) context. In this paper we apply it 
to the quantum many-body problem in TL, in particular 
in the Kicked Ising model (KI), namely the Ising spin 
1/2 chain periodically kicked with a tilted homogeneous 
magnetic field. KI is particularly interesting since it pos- 
sesses parameter-space regions with positive measure of 
non-ergodic behavior in TL surrounding the integrable 
cases JTo| of vanishing measure, which is an additional 
evidence for a conjecture on existence of intermedi- 
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ate, non-integrable and non-ergodic quantum motion of 
disorderless interacting many-body systems in TL. 

We start by rewriting the fidelity (Q) in terms of 
Heisenberg evolution of the perturbation A t := U~ t AU t 



t-i 



F(t) = (e 



iAq8 iA\8 



JA t -i§\ 



f{l[exp(iA t ,5)) (2) 

t'=0 



which is achieved by t insertions of the unity U U and 
recognizing l/ — (* "^UgU 1 = exp(iSAt'-x). T is a left- 
to-right time ordering. Next we make an expansion in S 
expressing the fidelity in terms of correlation functions 



F{ f) = l+ Y. l —r f Y. (A tl A t2 ---A tm ). (3) 

m—1 ti ,t2 .. .t m — 

Being interested mainly in the absolute value \F(t)\ we 
will in the following choose perturbations with vanish- 
ing first moment a :— (1/t) X)t^o(^*') — so that 
the series ([|) starts at m — 2, since a shift by a mul- 
tiple of unity A — > A — al simply rotates the fidelity 
F(t) — » cxp(— iaS)F(t). On the other hand, we can elim- 
inate not only the first, m = 1, but all odd orders in 
the expansion (0) by considering the symmetrized fidelity 
F{t) = {Ug^U s ^ 2 ). To second order in S we have 



(4) 



where it is assumed that 2-point time correlation function 
is homogeneous CU(t' — t) :— (A t A t i), as is the case for 
uniform average over initial states (.) = tr(.)/N. Eq. (g) 
reveals a simple general rule: the stronger correlation de- 
cay, the slower is decay in fidelity and vice versa. Below 
we discuss two different cases in the limit M — > oo : 

/. Ergodicity and fast mixing. Here we assume that 
CU(t) — > sufficiently fast that the total sum converges, 
S A ■= (1/2) Et=-oo C A (t), \S A \ < oo. For times t much 
larger than the so-called mixing time scale t 3> £ m ix which 
effectively characterizes the correlation decay, e.g. i m i X = 
Y, t \ tC A{t)\/ Et \ C A(t)\, it follows that the fidelity drops 
linearly in time F c (t) = 1 — t/r c + 0(S 3 ) on a scale 



t = S A 5 



(5) 



In order to show even stronger result we further as- 
sume fast mixing with respect to product observables 
B tt > = A t A t , with {B tt >) = C A (t' -t), of order k > 2, 
namely (B tlt2 B t3U ■ ■ ■ B t2k _ lt2k ) -> n^UC^-its,) a s 
ti,t2, . . • are ordered and £2.7+1 — — * oo. Therefore, 
the leading contribution for large t to each m-term of (||) 
comes from sequences (ii, t2, . . . t m ) where consecutive 

pairs (t2j-x, t2j) are close to each other, t^—tij-x ~ i m ix- 
Since for odd m time indices cannot be paired these terms 
should vanish asymptotically (as t — + oo) relatively to 
even m terms. Thus we can evaluate (2k — 1)!! equivalent 



sums over t'- = t 2 j — tzj-i giving, for t ^> t 



even m = 2k terms in Eq. (||) as fc-tuple of independent 
rms over t'- - 

oo 



{-l) k (2k-l)\\2 k 5 2k S k A 



k=0 



(2k)\ 



exp(-t/r e ). (6) 



Note that formulae remain valid in a more gen- 

eral case of inhomogeneous time correlations where one 
should take S A ■= Um t ^. 0O (l/i) Y,!°v=o{ B tt')- 

II. Non- ergodicity. Here we assume that auto- 
correlation function of the perturbation does not de- 
cay asymptotically but has a non- vanishing time-average, 
Da :— limt_,oo(l/t) Y?t~=o *7a(*')i though the first mo- 
ment is vanishing (A) — 0. For times t larger than 
the averaging time i avo in which a finite time average 
effectively relaxes into the stationary value Da, we can 
write fidelity to second order which decays quadratically 
in time, F nc (t) = 1 - (l/2)(t/r nc ) 2 + 0(S 2 ), on a scale 



t - D- l/2 c 

'ne — *-> A L 



(7) 



More general result can be formulated in terms of a time 
averaged operator A := limt_ >CJO (l/t) ^2*7— A t t, namely 
for t ^> i avc Eq. (H) can be rewritten as 



F nc (t) = l+J2 



ml 



■{A m ) = (exp(iA6t)). (8) 



Global behavior of F nc (t) for non-ergodic systems, where 
higher m-terms of (g) become important, depends gen- 
erally on the full sequence of moments (A m ). We argue 
below, by giving an example of spin 1/2 chains, that there 
are large classes of perturbing operators where these mo- 
ments can be shown to possess normal gaussian behav- 
ior, yielding Eq. (g). Non-ergodic behavior is certainly 
present for generic observables in completely integrable 
systems where a sequence of conservation laws can be 
used to estimate the time-averaged correlator Da [ pr{ , 
but we wish to make a stronger statement, namely that 
there is a generic regime of intermediate dynamics in non- 
integrable systems displaying non-ergodic behavior 

Let us now apply our theory to quantum spin- 1/2 
chains described by Pauli operators a xyz on a periodic 
lattice of size L, j + L = j, acting on a Hilbert space of 
dimension TV = 2 L , fix the average (.) = tr(.)/7V, and 
assume that our Floquet-operator U is translationally in- 
variant (TI) on a lattice. It is useful to introduce a set of 
local TI observables Z L = i" 1 / 2 J2j tf tf+i ' ' ' of 
order n < L, where s = [so, si . . . s n ], so, s n G {x, y, z}, 



G {0,x,y,z},l < j < n 



1, and er° 



1. Using 



( a j a k) ~ ^j.k^s,r one may derive a contraction formula 

U{a,/3}={1...2fe} 
all pairings a,0 

while for odd number (Z^Z^ ■ ■ ■ Z l2k+1 ) = C(L _1 ), 
hence Z s become independent gaussian field variables 
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in TL depending on a multi-index s of variable but fi- 
nite length. Therefore, any TI pseudo-local (PL) ob- 
servable A, having by definition Q Z 2 -expansion in 
the basis Zg_ (when L = oo), namely A = ^2 s as_Zs_, 
(A 2 ) — J2 S \ a ±\ 2 < 00 > possesses normal gaussian mo- 
ments (A 2f ) = {2k - l)!!(A 2 ) fc (l + 0{L- 1 )). Further, 
for a general TI PL observable A, its time average A is 
also TI PL, since it can be formally expanded in terms 
of Z s due to construction of A, and such expansion is 
I 2 since (A 2 ) = (AA) = D A < (A 2 ) @. However, 
for a more general non-TI PL observable A, i.e. such 
that its linear projection to the space of TI observables 

En=o ^loWoj+n * s P-^' one cann °t generally show 
that A is TI PL although we believe that this is a typi- 
cal situation, which we can prove in two cases: (i) If the 
spectrum of propagator U is non-degenerate (for any fi- 
nite L), then the matrix of A is diagonal in the eigenbasis 
of U and A is TI due to Bloch theorem, (ii) If the system 
is integrable having a complete set of TI PL conservation 
laws Q n , n = 1, 2 ... in the sense that {Q n } is a complete 
set of eigenvectors of the Heisenberg map IAA = U'AU 
for eigenvalue 1 then the time average is a projection 
A = J2n(Qn A )Qn (assuming that (Q n Q m ) = 8 nm ) 
which is TI PL. This is the case for KI model studied 
below. Finally, assuming either (i), (ii), or simply TI PL 
perturbation A, we find that moments of time-average A 
are gaussian (A 2k ) = (2k - l)WD k A {\ + C(L -1 )). Sum- 
ming up the formula (ra) produces gaussian decay 

F nc (t) = exp(-(i/r„e) 2 /2), (9) 

for t ^> i avc , on a time scale (0), which can be computed 
in a typical integrable situation (ii) as shown bellow. 

Few remarks on the case of finite dimension Af < 00 
are in order: (1) F(t) will then start fluctuating around 
zero with magnitude Ffl uct = AA -1 / 2 for very long times 
t > t* (Af) where the time scale t* (Af) is determined from 
the condition F(t*)y =oc = Af~ 1/2 . (2) F(t) decays all 
the way down to TV -1 / 2 only for a typical or random ini- 
tial state \ip) with ~ Af non- vanishing components when 
expanded in the eigenbasis of U, or for an average over 
\ip). If on the other hand one considers the initial state 
which, when expanded either in the eigenbasis of U or of 
Us, contains essentially only few, say m dominating com- 
ponents, like the regular coherent state of Peres ||, then 
F(t) is a quasi-periodic function with m small frequencies 
oc 6 and amplitudes ~ 1/m. (3) Even in asymptotically 
ergodic situation the correlation CU(t) has a plateau for 
finite Af, which can be estimated using random matrix 
model for the propagator U % as Da ~ D* A (N) := caIM 
where ca is some constant with respect to M . The non- 
vanishing correlation plateau gives a dominant contri- 
bution to Eq. (Q) resulting in a quadratic (or gaussian) 
decay of F{t) as soon as t c > Sa\m=oo/D* a , i.e. when 

5 < Sp(J\f) := S~^~ c~J^ N~ x l 2 . This perturbative regime 
of very small perturbation strength, existing for finite Af 
only, is consistent with the first order perturbation ex- 



pansion of eigenstates of Us in terms of the eigenbasis of 
U @. 
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FIG. 1. Correlation decay for three cases of KI: (a) in- 
tegrable h z — 0, (b) intermediate h z — 0.4, and (c) er- 
godic h z = 1.4, for different sizes L = 20, 16, 12 (solid-dotted 
connected curves, almost indistinguishable in (a,b)). Circles 
(a) show exact L — 00 result. Chain lines are theoreti- 
cal/suggested asymptotics (see text). 

Consider an example of KI model with the hamiltonian 
L-i 

H K i(t) - i J *°j°j+i + S P (t)(h^j + h,a*)} (10) 
3=0 

where S p (t) = J2 m — m P)i with a Floquet-map U — 
exp ( - i J z <r| cr| + 1 ) exp ( - i Y,j ■ ( K aj + h z a?)), where 
we take units such that p = h = 1, depending on a 
triple of independent parameters (J Zl h x ,h z ). KI is inte- 
grable for longitudinal (h x = 0) and transverse (h z = 0) 
fields [0, and has finite parameter regions of ergodic 
and non-ergodic behaviors for a tilted field (see Fig. 1). 
The non-trivial integrability of a transverse kicking field, 
which somehow inherits the solvable dynamics of its well- 
known autonomous version Jlif , is quite remarkable since 
it was shown jl0| that the Heisenberg dynamics can be 
calculated explicitly for observables which are bilinear in 

fermi operators Cj = — icr|) 11/ <J °/ with time cor- 
relations decaying to the non-ergodic stationary values 
as \C A (t) - D A \ ~ t- 3 ' 2 Q. For D A we find explicit 
expressions, the simplest, 
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D„ 



max{| cos(2J z )|, | cos(2h x )\} — cos 2 (2h x ) 
sm 2 (2h x ) 



(11) 



and Dm = LD a *, for the component of spin crj, and the 
component of magnetization M = J^j a j > respectively. 

In a general situation of non-integrable KI we wish to 
test our theory by a numerical experiment. We consider 
a line in 3d parameter space with fixed J = l,h x = 1.4 
and varying h z exhibiting all different types of dynam- 
ics: (a) h z = integrable, (b) h z — 0.4 intermedi- 
ate (non-integrable and non-ergodic), and (c) h z = 1.4 
ergodic and mixing. In all cases we fix the operator 
A = M which generates the perturbation of KI model 
with h x —> h x + (hi + h 2 z hcoth)S/h 2 + 0(S 2 ), h z -> 
h z + h x h z (l-h cot h)S/h 2 + 0(S 2 ), where h = yfh 2 +h 2 , 
and vary L and S. Since we want the perturbation 
strength to be size L-independent we scale it by fix- 
ing 5' = S^L/Lq where Lo := 24. Time evolution 
has been computed efficiently by iterating the factored 
Floquet map (in terms of 1-spin and 2-spin propagators 
- 'quantum gates'), requiring oc L2 L computer opera- 
tions per iteration per initial state. In integrable case (a) 
we confirm saturation of correlations to the theoretical 
value (To) D M = 0.485126 x L (Fig. la), as well as gaus- 
sian decay of fidelity (^|) with time-scale r ne given by (0) 
which terminates at t ~ t* c = -^(lnA/") 1 / 2 (Fig. 2a) In 
non-integrable (intermediate) case (b) , we find persisting 
non-ergodic and non-mixing behavior since rescaled cor- 
relation functions of typical observables CU(£)/ (A 2 ) relax 
on a short L-independent time scale to a non-vanishing 
value Da/ (A 2 ) and converge to TL very quickly with 
increasing size L (Fig. lb), but as opposed to inte- 
grable case (a) the relaxation appears to be exponential 
\C M (t) - D M \/L ~ exp(-i/i avo ) with t avc w 7.2 (inset 
lb). Such behavior has been observed for other two com- 
ponents of the magnetization M V ,M Z and supports ex- 
istence of intermediate dynamics observed previously in 
kicked t-V model |]. In Fig. 2b we confirm gaussian de- 
cay of F(t) predicted (0) from numerically observed value 
of Dm — 0.293 x L, again up to time i* e (2 L ). In ergodic 
case (c) we find fast decay of correlation functions fitting 
well to an exponential \CM{t)\/L ~ exp(— t/t m i x ), with 
imix ~ 6.0. Consequently we find exponential decay of 
F(t) of eqs. (§,§) using S M = (1/2) £ t C M (t) « 2.54 xL, 
up to the saturation time t* = (l/2)r G mTV (Fig. 2c). 

In conclusion, we have presented a simple theory for 
the stability of quantum motion with respect to a static 
perturbation of the evolution operator in the limit of 
Hilbert space dimension TV — > oo, characterized by the 
fidelity measuring the distance between time evolving 
states. The fidelity was expressed in terms of inte- 
grated time-correlation functions of the perturbing op- 
erator, showing that faster decay of correlations gives 
slower decay of fidelity, meaning that 'chaotic' dynamics 
is more stable in Hilbert space than 'regular' one (unless 
the state that one is looking at is simply related to the 
eigenstates of the system)! In the two limiting cases of 
mixing and integrable (or more generally, non-ergodic) 



dynamics we find, respectively, exponential and gaussian 
decay. For example, our finding has strong implication 
for the stability of quantum computation with respect 
to static imperfections (e.g. uncontrolable residual in- 
teraction among qubits) ]l5] . In other words, Eq. (j^) 
is a version of the fluctuation- dissipation formula for the 
'dissipation coefficient' l/r c of Eq. (JsJ) which diverges in 
non-ergodic regime. If the system has a well defined clas- 
sical limit then our formula (|^) has a clear and simple 
classical limit h — > too, with an integrated classical au- 
tocorrelation function substituting the quantum one p6| . 
We speculate that our finding is a manifestation of "the 
structural invariance" |17| of quantum chaotic dynam- 
ics. Although in this paper our theory has been demon- 
strated in a specific kicked many-body problem, namely 
the quantum kicked Ising spin 1/2 chain, we should em- 
phasize that it should be generally valid (within the time 
and perturbation scales depending on the Hilbert space 
dimension) and thus applicable to any unitary evolution, 
in particular also to any experimentally interesting quan- 
tum dynamics. 
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FIG. 2. Absolute fidelity \F(t)\ for three cases of KI: (a) 
integrable h z = 0, (b) intermediate h z = 0.4, and (c) ergodic 
h z = 1.4, for different sizes L — 20, 16, 12 and different scaled 
perturbations 8' . Chain curves give theoretical predictions. 
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